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ABSTRACT 


A series of Hieasurements of near field pressures and. 
turbulent- velocity fluctuations were made in a low speed jet 
with a Reynolds number near 50^000 in order to investigate 
more quantitatively the character and behavior of the large 
scale structures and their interactions with each other. 

Time averaging techniques and simple eduction 
techniques were insufficient to detect interactions betv/een 
vortex structures beyond a diameter dovmstream. . More 
sophisticated detection schemes were investigated with 
little success. Finally, the near field laeasurements were 
modelled according to the "vortex pairing” hypothesis to' 
deduce the distribution of pairings along the jet axis and 

■■ , t- 

the variances about the mean locations . 

The hodograph plane description of turbulence vras 
explored in some detail , and a complex correlation quantity’ 
was synthesized which has useful properties for turbulence- 
in the presence of mean shear. | i 
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1 . INTRODUCTION 

! 

\ 

1.1 Motivation 

There is substantial evidence that there are large 

scale organized structures within the jet mixing layer. 
They are large scale in the sense that they are spatially 
coherent over a major fraction of the mixing layer 

thickness, and they are organized in the sense that their 

associated flow field is induced by inherently compact 
regions of concentrated vorticity. We feel that their 

existence has been adequately demonstrated at least during 
the initial stages of shear layer growth. V!e are interested 
in the contribution of these structures to turbulent mixing 
beyond the potential core as well as their’ contribution to 
the radiated noise. 

In the past, the interpretation has been that noise is 
radiated by local fluctuations of momentum flux, and the’ 
resultant acoustic radiation is the linear superposition of 
these purely local contr ibutions . It would be of great 

interest to examine whether these fluctuating momentum' 
fluxes are associated primarily v;ith interactions between^ 
the large structures . There are a number of observations 
that seem to suggest this possibility: 


1 



2 


1) From flow visualization experiraents conducted by' 
Browand, Chu, and Laufer (1975), it is apparent that at 
least for moderate Reynolds numbers (Re < 10,000), the 
initial growth of the jet mixing layer is dominated by 
vortex ring structures that coalesce as they interact. The- 
coalescence is sudden and results in an explosive ©reduction 
of smaller scale radial and azimuthal flow. _ Both thJ 

exchange of momentum between rings and the radial 

* 

accelerations are plausible noise machanisnis. ; 

2) A collection of total acoustic power spectra] 
compiled by Howes (1960) and discussed in detail below, 
indicates that for large Reynolds numbers and moderate Mach 
numbers (Re > 100,000 and M > .6) the spectral peak scales 

with jet diameter independent of viscous, compressible, or! 

* 

thermal influences. This scaling is consistent with the 
observed ring structures since their size and spacing is 
intrinsically related to jet geometry. ^ 

• , ’ i ■ 

3) Recent shear layer measurements of Winant and 

Browand (1974) indicate large intermittent increases in 

■ ' ' ' ' ' , , ' ' • ■ ' ' ■ ■ ■ ■ \ ■ 

momentum flux associated with passage of vertex structures 

and wi th interaction of adj acent structures . Intermittent 
momentum fluxes, after all , are the contributing noise 
sources according to the Lighthill (1952) formulation. 







4) Powell (1964) has shown that when a free vortex 
ring stretches it radiates acoustically like a dipole/ with 
a far field amplitude proportional to the second time 
derivative of its Kelvin impulse. The kinematics of vortex 
coalescence suggests that impulse is exchanged between 
rings. In that case vortex coalescence would appear 
acoustically like a longitudinal quadrupole radiator. 

though' this evidence is at best indicative/ it 
suggests a conclusive experiment: if it were possible to 

detect vortex coalescence v/ith probes in the shear layer or 
near field of the jet/ then simultaneous directional far^ 
field measurements V 70 uld establish any causal relationship.; 
Detecting vortex coalescence in a high Reynolds number^ 
compressible jet using a small number of probes is a 
formidable task. Consequently/ this work is an 

investigation of techniques aimed at detecting such events. , 

I 

1.2 Parametric survey of jet noxse spectra 

■- i- 

For the present investigation a low speed (100 fps) 
moderate Reynolds number (50,000) jet was selected. This^ 
jet seemed a reasonable compromise xn that flow, 
visualization experiments have demonstrated vortex pairing’ 
up to a Reynolds number of about 10/000, while noise 
measurements have been made on jets dov;n to a Reynolds 



nuiiiber of about 100,000. Since the jet is essentially 
incompressible at these low speeds, hot wires respond 
directly to flow velocity. 

The question remains whether measureiaents made on a low 
Reynolds number, incompressible jet are applicable to a high 
speed, compressible jet. A \9ay of approaching this question 
is to examine jet noise scaling with Reynolds number and 
Mach number. Such a parametric survey of total acoustic 
power spectra, derived from a compilation by Howes (I960) 
are summarized in figures 1 through 4. The spectra 
represent a v;ide rang?:, of geometries and experimental 
conditions. They include both air jets and engines, 
exhausting through both parallel and convergent nozzles. 
The test facilities included reverberant chambers 
(Fitzpatrick and Lee, 1952; Tyler and Perry, 1954; 
Waterhouse and Berendt, 1958), as well as free field 
surroundings (Lee, 1953; Rollin, 1958; and Howes, 1960). 

Figures 1 and 2 show acoustic power spectra as a 
function of jet Strouhal number (fD/Ue), with Reynolds 
number (UgD/v) as the parameter, for both moderate and high 
Mach nximber jets (M = . 60 and . 90) . Exit Mach number , based 
on local sonic velocity, is the parameter in figures 3 and 4 
for both moderate and high Reynolds number jets (Re === .20 
and 2.0 E 6) . The power spectra are normalized to unity 
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Figure 1. Parametric study of jet noise spectra 
M = .6 and Re = .2-2.0x10 
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Figure 2. Parametric study of jet no i s 
K - .9 and Re = .09“3*3xl0^ 

















area, so the relative levels are less importaht than the 
locations of the spectral peaks, and the slopes of the high^ 
and low frequency asymptotes. The effect of Reynolds number- 
on spectral similarity is clear: there is a consistent^ 

viscous influence in that the high frequency asymptotes' _ 
flatten with increasing Reynolds number while the low^ 
frequency asymptotes steepen, implying a transfer of 
acoustic energy from larger to shorter wavelengths. The 

value of the spectral peak , however , remains relatively^ 
constant with Reynolds number. If the spectral peak 
represents the characteristic time scale of the radiating; 
get structure, the implication is that the structure scales; 
independently of Reynolds nuiaber. Figures 3 and 4 indicate^ 
that the effect of compressibility on acoustic scaling is^ 
less certain. There does appear to be some shift towards 

smaller scales at the lowest Hach numbers. This may simply 
indicate a departure from purely quadrupole radiation, for| 
instance by acoustic sources interior to the nozzle as has^ 
been often suggested for explaining departures from the^ 

eight power law at low velocities. ^ j 

■ ■ , ■ ■ ■ ■ ■ } 

■ i 

. .. . . . ■ ■ . .. ■ ■■■ ■ •■...■■ ■ ■ .1 

.1,3 Review of flow visual izarion s uudies | 

The existence of vortex rings and their mutual 

interaction in a low Reynolds number Cl 0 ,000) jet were 

observed as early as 1955 by k. 3. C.^ 



interested in the physics of pipe tones produced in a jet 
issuing through an orifice plate, and use^ a double spark* 
shadowgraph to visualize the jet structure. He observed 
that the primary vortex shedding frequency produced the 
fundamental pipe tone, and the pairing mechanism excited 
various subhannonics. 

Rockwell and Niccolls ( 1972) recognized \-ortex pairing 
as the mechanism for the shear layer growth in the potential 
core region of a two-dimensional jet. They used the 
hydrogen bubble technique to visualize the flow, and 

included a Reynolds number range of 1860 to 10,800. Both 
symmetrical and asymmetrical pairings were observed. 

That vortex pairing is the mechanism responsible for 
layer growth was amplified . by Winant and Browand 
(1974) in their flow visualization work vrith a 

two-dimensional, turbulent mixing layer. They injected dye 
into the mixing layer formed between two parallel streams of 
water moving at different speeds. In particular they 
observed that the onset of pairing marked the tr^nsif ion of 
the shear layer from laminar to turbulent grov;\-h arid, that 
■yortex pairing inherently produces linear shear layer 
spreading . 

Similar dye studies by B‘fov;and, Chu, and Laufer (1975) 
in a circular water jet demonstrated that the spreading rate 



of the mixing layer can be increased by artificially driving 
the jet core at a subharmonic of its initial instability. 
Forcing seems to increase the spatial coherence of the flow 
and perhaps hastens the transition from laminar to turbulenw 
spreading . This study v;as conducted at j et Reynolds numbers 
up to 15,000. 

1.4 Review of mixing layer measurements 

Flow visualization studies have proved . invaluabre in 
providing a physical model but ate fundamentally qualiuauive 
techniques, and are confined to low and moderate Reynolds 
number flows. Quantitative turbulence measurements in 
moderate to high Reynolds nximber jets have until recently 
lacked a direct connection with the flow visualization work. 

j. C. Laurence .(1955) compiled a statistical survey 
of the mixing layer of a 3 . 5 -inch circular jet. He used 
single wire and x-wire probes to measure .moments and 
correlations at varying radial positions along the potential 
core and transitional mixing regions of the jet. The 
correlation measurements were used to dea.ine lateral and 
longitudinal length scales, which were observed to be 
insensitive to Reynolds number over a range 37,500 to 
350,000, and insensitive to Mach nxmber over a range .2 to 
7, He found that the peak turbulence level occurs at the 


radial position r/D =' . 5 . ! 

Hot V 7 ire measurements by Bradshaw, Ferris, and Johnson 
( 1964 ) indicated that irrotational velocity i.lucuuations 
outside the jet near the nozzle are due to large scale, 
eddies within tlie mixing layer. From measurements of the 
correlation v^w^> A v^><w^> they dete3rmined rhat the radial^ 
(v) and azimuthal (w) components are generally either in 
phase or in antiphase, implying that the flov7 assocxated 
with the eddies is not axisymmetric, but that the large 
eddies are characterized by random swirl evenly distributed 
so that the mean <v;v> - 0. ' . ; 

Mollo-Christensen ( 1967 ) was interested in the. 
contribution of the large eddies to the near field static. 
pressure fluctuations. He placed microphones at the edge of 
a . low speed jet and made both broadband and narrowband 

crosscorrelations between microphones that were separated, 
axially and azimuthal ly along the potential core. From 
these measurements he postulated an intermittent, spatially 
regular structure with random phase. Furthermore the 

frequency of the regular structure was observed to shift 

discontinuous ly to its own subharmonic -with downstream 
distance. From comparisons between crosscorrelations of 
narrowband pressure signals %^hich had been filtered about 
different center frequencies, it was in^errea uhat the 
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transfer of energy v/ithin a wave packet from smaller to 
larger wavelengths was a nonlinear process. I 

. Mollo-Christensen' s observations were supported by Ko 

• I . ! 

and Davies (1 971) who made hot wrre measurements inside the 
potential core of a low speed jet over a velocity range 
tJe = 165 “• 335 fps. They found that the velocity! 
fluctuations in the core vrere characterized by a relatively 

1 ' ^ i 

narrow power spectrum and a mean convection speed ratio 
;Uc/Ue = .60 - .65, independent of exit speed Ug. The: 

i ... ■ . ■ I 

% 

|convection speed is so low that the velocity fluctuations 
icould not be carried by the local mean flow, implying that 
the fluctuations are induced from v^ithin the mixing layer.’ 
'As a spatially integrated effect the induced fluctuations 

I ^ i 

would be most sensitive to organized structure even if the 
latter did not dominate the turbulence locally. ■ 

. Crow and Champagne (1971) investigated the influence of 
initial conditions on the mixing- layer structure over a jet 
Reynolds number range 62,000 to 124,000. Their jet plenum 
chamber v;as driven acoustically at various harmonic 
resonances’, introducing fluctuations into the core v/ith 
amplitudes as large as 555 of the exit velocity Ug . When the 
jet was driven at a Strouhal number of .3, which had been 
observed to be the dominant node, the resulting disturbance, 
V7hen filtered around St = .3, grew exponentially in the 



‘first couple of diameters downstream, reacned a peak level 
Uar the end of the potential core, and thereafter declined, 
icrow and Champagne interpreted this in terms of a linear 
Iwave instability followed by non-linear saturation. However 
Uen they drove the jet at a Strouhal number of .6, the - 3 
istrouhai disturbance was driven more effioxenrly than at ius 
•own wave number. This apparent subharmonxc amplxfxcatxon is 
iconsistent with vortex pairing , which can also account for 

■the observed non-linear saturation. Crow and Champagne 

I 

ifound the jet shear layer to be non-dispersive at 


■jfrequencies greater than Strouhal .3, the dominant mode, and 
^dispersive at smaller wavenxambers , raising a possible 
connection between disper s ion and the observea concentration 
'of turbulent energy about a dominant Strouhal number. 


- Lau, Fisher, and Fuchs (1972) identified the mean 
structure within a circular jet shear layer at axial 
position x/D=2 as a series of vortex rings moving at a mean 
^speed U^/Ug =• .6, and with a mean spacing of 1.2 diameters. 
They used x-wire probes to measure velocity fluctuations 
inside the core, in the entrainment region, and inside t_he 
mixing layer of their jet which v;as operated at a Reynolds 
‘number of 210,000. The vortex model v;as deduced from phase 
reiationships between radial and axial velocity rluctuations 
measured in the entrainment and core regions. Correlation 
measurements _betw-een a. microphone, and the hot wire „j>robes 


'indicated that vortex passage manifests itself in the near 
field as a static pressure trough. In particular they 

I 

proved for the first time that the near field pressure 

jfluctuations are produced by the passage of vortex rings 

j ■ ■ ■ , , ^ ■ ■ 

'■through the mixing layer. 

I . ■ . . . 

! Davies et al (1974) were particularly interested in 

obtaining quantitative information about the interactions 
•between vortex rings that are so evident in the flow 

Ivisualization studies . Arrays of single and x-wire probes 

.1 ■ ' ■ 

Vere placed inside the mixing layer and space-time 
^correlation as v;ell as eduction techniques vrere utilized to 
bring out the interactions. At moderate Reynolds numbers 
i(Re = 20,000) they were successful in detecting one or two 
'coalescences near the nozzle exit, 'which they called 
!•Marainar" pairings. However vortex pairings could not be 

i 

'detected further downstream than a couple of diameters, and 

i .... 

-. I ■ ■ ■■ ■ ■ ■ ■ ■ ■■■ 

'even the laminar pairings could not be detected at higher 
jReynolds numbers (Re = 200,000) •. From flow visualization 
^studies of their own, using smoke jets, and individual smoke 

I"; .' ■ ' ■ 

rings, strong azimuthal structure was demonstrated. In 
particular there vtere occasional events that appeared like 
secondary, radial jets that are possible noise mechanisms. 



'l.5 Scope of present work 

j This study is an investigation of techniques aimed at 

detecting organized, intermittent events inside the mixing 
* ■ " , , 1 
layer of a moderate Reynolds nxiraber circular jet using a 

ismall number of probes. Particular emphasis is placed on 

pressure signatures in the near field, and on the use of 

■ 1 ■ 

'conditional statistics to describe the mixing layer 

iturbulence. 

i ' ■ 

I The investigation was only moderately successful in 

detecting interactions betv/een vortex rings, and v;as 

^unsuccessful in detecting any organized structure at all in 

I ' ■ ^ 

;the transitional mixing region. However, based on the 
measurements, an indirect case can be made suppbrting the 

i ■ . ■ 

Ivortex pairing hypothesis as the mechanism responsible for 
'the linear spreading of the jet shear layer just as it is 
ithe mechanism for free mixing layer growth. Chapter 7 is a 

I 

isummary of this evidence, while overall conclusions appear 
'in chapter 8. The actual near field and mixing layer 
jmeasurements are discussed separately in chapters 5 and 6/ 
irespectively. The special statistical techniques used in 
•obtaining these rrieasureinents are discussed in chapter 3. 
Chapter 4 contains a preliminary survey of the jet in terms 
of profiles of mean velocity and rms turbulence, which is 
intended to establish the model jet as normal. The problems 



17 


:and procedures involved in calibrating the x~wire probes are 
discussed in the appendix. 




2. FACILITIES 


2.1 Model jet 

I 

I 

Most of the experiments described in this investigation, 
have : been carried out in the low speed modei. jet shown in 

■ I ■ ' ■ • 

figure 5 . It consists of a vacuum cleaner blower , a 

settling chamber xjfith dsimping screens, ana a circular 
nozzle. The contraction ratio of the nozzle is 36:1 and the 
exit diameter is 1.00 inch. The maximim attainable exit 
velocity is 220 fps, but most of the work was carried out at 
velocities near 100 fps, corresponding to a Reynolds number 
(UeD/v) of 52,000. 

The exit boundary layer is laminar and xts thickness is 
measured to be .05 diameter. The exit plane mean velocity 
profile is flat to better than 0.5^ and the exit plane 
turbulence level is 0.3^ of exit velocity Ug. This 
relatively high level was considered satisfactory for the 
type of measurements planned. 

2 . 2 Instrumentation 

Mean velocity measurements v;ere made wxth a total 
pressure Pitot tube. Depending on the application, total 
pressure V 7 as either read xrom a BETZ manomeuer, or recorded 

18 '■ i ■ • 
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off of a VALIDYNB Model DP15TL capacitance type differential 
pressure transducer. 

t 

Near field pressure measurements vrere made using 
BRUEL 6 KJAER model 4138 eighth— inch, capacitor microphones 
with a nominal sensitivities of -SO dB re 1 voit/ubar, and a 

i 

nominal bandwidth of 20 Hz to 150 kHz at levels below 170 

I . 

■ r ^ 

dB. Microphone output was amplified through, a BRUEL S KJAER 
model 2618 cathode follower. 

: I * 1 ^ ^ ^ ^ 

i . - ■ ' • ’ ■ 

] Turbulence maasurements v/ere made with constant 

i ^ 

pemparature hot wire anemometers. The wires were .0001 inch 

diameter Platinum, with an aspect ratio of about 400. The 

* - . - 

bridge circuit, shown in figure 6 , v^as designed for an 
overheat ratio of 2 . In the: presence of a large ratio • of 

mean velocity to r.m.s. fluctuating velocity, the wire 

! ' ; ■ ’ ■ 

responds to the velocity component parallel to the mean 

' I ... 

fiov;. Consequently, a single wire is sufficient to record 
! • • . . . ' ' 
axial velocity fluctuations. However, an x~wire probe is 

needed to decompose the turbulence into axial and radial 

components. It is necessary to determine calibration 

polynomials in order to obtain orthogonal velocity 

components from the outputs of the two wires . As particular 

care must be exercised in deriving and measuring these 

polynomials, the calibration procedure is described in 

detail in the appendix. . 
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A traversing machanisra of our own design was used for 
radial positioning of probes, and the model jet itself was 
moved for axial positioning, i'Then greater precision V7as 

required, a PRIOR three-dimensional laicropositioner was 

t ■ ■ ■ _ . ■ 

used; with this instrument a probe can be positioned to 
within .1 mm. ; 

I ^ i 

i ■ ; 

|2.3 Digital equipment • 

’ ... \ 

I Most of the correlations and educted signals v7ere 

i ^ - : 

measured with a HEWLETT-? ACKAPvD model 3721 A digital 

correlator . This device has the capahilities of generating 

'autocorrelations, crosscorrelations, and histographs , as 

well as signal recovery provided an appropriate trigger can 

.be generated. The correlator has a sampling rate that is 

adjustable in steps from 1 Hz to 1 MHz . It was used for 

most of the near field pressure measurements and for the 

single v/ire turbulence measurements . 

; For measurements of joint statistics , we chose to 

•linearize the x-wire signals digitally. The raw signals 
were digitized using a subsystem, of our own design, 
described by Jentes (1975) . It consists of 16 sample and 
hold circuits and an analog to digital converter interfaced 
with an INTERDATA 5 programmable lainicomputer . The maximum 
sampling rate in microseconds is 100 + 20* (= channels) . the 



digitized record is vrritten onLo disk, then transferred to 
digital tape. 

A PDP 10 general purpose computer was used for most of 
the data reduction and analysis. 


; 3. STATISTICAL TECHNIQUES 

3.1 Correlation techniques 

I Correlation aeasureinents have been the traditional 

means of determining time and length scales inside mixing 

! •: ■ 

Layers . The tvra point space-time correlation is defined: 

1 ■ ■ 

' ^ ■ J 

I <e,(x)e„{v)>/ = liin_JLj e,(x,t)e^(y,t+T) dt (1) 

i 1 - Z - u; - 21 -T 

' ■ 

and in general contains information about the spectral 
icon tent of the turbulence, the coherence length scales, and 

t ■ 

_ t , . ■ 

jthe mean convection speed of the turbulence . 

! ■ . ' . . 

I The coordinate system is defined in figure 7 , and will 

■f - ■ ■ 

be referred to throughout this paper* 

i ■ ■ 

. i ■ , , , . , ■ ■ ; ■ ' * 

\ , ■ ■ ■■ ■ "•■■■■■ ■ 

I 

3.2 Two point eduction 

,i ■ ' ■ , ■ ■ ■ , ■ ■ ■ ■ , ' ■ ■ '' '■ 

• Two point eduction is a type of flow visualization, and 

'^Qj^glgts of generatxng a trigger pulse from conditions at 

: i ; i, ; : ; : • ;i _ / . 

lone point in the flow field and sampling the flow at' a 
'second location after some fixed time delay relative to the 
pulse- When the sampling is done over an ensemble, the 
■average as a function of time delay is referred to as the 
"educted" signal; a number of educted signals at various 
positions in the flow field determine the relationship 

' 2 k 
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between any spatial structures in the field and the 
triggering condition at the triggering location. 


I The most difficult aspect of tvra point eduction is in 
choosing an appropriate triggering condition. There are two 
main considerations in selecting the condition: 
i 1) it should be exclusively characteristic of the 

; flow field structure one is tryi.ng to visualize, 

; 2) and it must be temporally narrow to provide a 

t' ’ ■ ^ \ ^ ^ 


i precise time reference . 

The trigger must be derived from some feature of the basic 

\ 

random signals. The simplest form of such a criterion is 

! ’ i 

-Signal slope and level, . in much the same manner as an 

bscilloscope trigger signal is formed. Kowever such an 

jinclusive condition can bring out only very regular 

structures - In such cases eduction provides little more 


Information than conventional crosscorrelations . The real 
advantage in eduction techniques lies in detecting 
subprocesses' , such as interactions between regular 

t ■ ■ , , . '/ ■ " ’ , 

structures , that are . intermittent, and lack a consistent 
phase relationship to the mean structure . Obviously , the ’ 

better the subprocesses are understood, the easier it is to 
choose an exclusive trigger condition. • 




3.3 Joint statistics 

One potehtially useful method of generating a trigger 
condition/ borrowed from pattern recognition, is to study 
statistical properties of the flow at a point as it appears 
in the hodograph plane, with the intent of identifying 
subprocesses. The conceptual basis for such an approacn 
derives from the nature of the dynamical equations. If the 
flow is incompressible, the pressure term in the equation of 
motion can be expressed as a volume integral over the o.lov7 

I '■ ■ . . ■ ■ 

’field at that instant. Consequently, the field equation 
involves only first differences in time, which implies that 
Ithe entire field at a given time is completely described by 
ithe entire field at the previous instant in time. However, 

f 

•if the field is known only at a single point in space, it is 
’uncertain how the field at that point vzill change over the 
hext instant. It is plausible, however, that the sta t idi es 
of the field at a point depends only on the statistics of 
the field at that point in the previous instant. Gur 
approach has been to measure various first differencess pa- 
the field at a point, and to attempt to associate them with 
deterministic subprocesses. 

The two joint statistical huentities of interest are 
the joint probability density P(u,v) and the joint 

conditional probability density a The^^^^^^ 
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probability density/ vzhen suitably normalized, is the 

t. 

frciction of tiitie tins field occupies position iU/V) in the 
hodograph plana. Measured ?(u,v) is incetendant of the 
sampling rate provided the sampling is faster t_no.n i_he power 
'spectral bandwidth of the turbulence. The joint conditional 
probability density represents the probability that the 
jfield would occupy position (u,v) after tine step "dt". 
■provided it started from position (u,v) . The dependence of 
lQ(u,v;u,v) on sampling rate v/ill be discussed later in this 

section. • 

i- . . ■ ' . . ■ j 

■ i ■ , . r .. . 

' ' 

i The proper normalizations for the probability densities 

!' ■■ ■ ■ ■ ■ I 

■i . . ' . ■ . . ' 

jare * ■. ■ 

t . . 

l cb - - \ I 

i // P (u, v) dudv = 1 (2a) 

I ■ ■■ ■ — ' . ... 

t . — oo . ■ .... 

. I ■ . . 

•and; , . ' : • 

1 ■ ■ 

i CO 

j I! Q(u,v;u,v) dudv = 1 (2b) 

i .—.03 

t 

I. . 

5 • . ■ ■ • 

Time averages of the turbulence can be computed as moments 

of the joint probability distribution. Specifically, mean 

Velocities are defined ■ 

■ CO 

‘ . < g.> = //q- ?(u,v) dudv (3) 

. * I. ■ ■. ' .^C3. ■ . ■ ■ ■ 

variances are defined . 


aj2= /j{ q.2_<qj>2} P(u,v) dudv 

— CO 

and stresses are defined 

oo 

— 00 

The correlation coefficient is defined 

i - ■ 

; R = {< uV>-< U>^'v>}/u'v' 

t ■ ■ ■ , ■ 

Adhere q ' , the rms fluctuating component , is 
; • <Ji' 

I 
t 

I ■ . 

^ The hodograph plane description of turbulence is not a 

new technique. It was used for example by Liepmann and 
Laufer ( 1 947 ) . However , joint probability densities are 
difficult measurements using analog equipns^it/ and it hasn «_ 
been until the advent of the high, speed digital computer 
that these kinds of measurements have become practical . - In 
practice they are among the easiest calculations to program. 
Pre-recorded recorded x-wire signals e^ (t) and e 2 (t) are 
converted to parallel digital records using sample and hold 
circuits. The records are linearized using an appropriate 
calibration scheme, resulting in parallel records of axial 
(u) and radial (v) velocities. The hodograph plans is 
divided into the desired number of cells lin this case 2500 
( 2 giis) _..o.£- - dimension — ^u x ^ v , . and each occurrence — O j — ( u , v.) 


(h) 


(5) 


( 6 ) 


(7) 



wi-thin the range Vj^+^av/2) generates a unique index 

associated with cell (u^,Vn) . By counting the number of 
such occurrences a joint histogram is generated, which is 
normalized according to equation (2a) to produce tne joint 
probability density. 

At the same time the joint proba.biliuy is being 
generated it is easy to look forward one time step ”dt" to 
measure the average forward differen.ces du> , < dv> , < du > , 

< dv^> , < dudv> associated with each cell (u,v) . These 

average differences describe the mean properties of the 
joint conditional probabilities Q(u,v;u,v) . For example the 
vector- (u-K du> , v-j< dv> ) defines the xaean conditional 
(< u> ,< v> ) associated v;ith the cell (u,v} . -?he tensor Tjj 
defined by 

■ Tjj Cu,v) = (<dq.dqj> “ < <^^ 1 ^ , v) 

can be diagonalized to describe the principal axis rotation 
and second moments of each Q(u,V/U,v) . Tn^o uoe Oj. forward 
<31ff©rences to estimate Q is ill'ustrated in _igure 8, and 
saves the prohibitive amount of computer storage that would 
be necessary to directly measure the conditional probability 
associated v/ith each cell. Note that the magnitudes of 
first and principal second moments depend directly on the 
magnitude of the time step dt? hovJever if d— x=> su^j-iciently 
small, the direction of the first moment anc. tne principal 
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axis rotation of Tjj are invariant with saraplxng rate. 

From figure 8, the vector (<du>,<dv>) can be 

visualized as describing the 'drift' of the mean of the 
conditional probability density, while the x^ensor T.j(u,v) 
describes its 'diffusion'. If the magnitude of the drift is 
much larger than the eigenvalues of the diffusion, the 
motion associated with that position in the hodograph plane 
can be considered largely deterministic; if it is the other 
v/ay around the motion can be considered noise driven. Based 
oh this concept, one can formulate the time averaged drift: 

03 

DRIFT = // hq[^x hq|^> P(u,v> oudv (9) 

— Op 

and diffusion 

CO 

DIFFQSIOISI = // '/Tkic dudv (lO) 

— OP 

Then the ratio DIFFUSION/DRIFT is a measure of the overall 
■•noisiness' of the measured velocity field. 


3.4 Complex correlation 

There are two general properties of the drift vectors 
dq(u,v) that are predictable. First, the drift will be 
largest in magnitude from cells (u ,v) farthest from the mean 
(^ n> ,< v> ) . Second, they will generally point radially 
towards the probability peak, in the sense that the most 

QEIGMAIi 
^ POOR QtJAlffl 


33 


< 33 Txt*t wxix t)3 tow3XC3 th.© most st©t© • 

However, if there is some intermittent subprocess present 
with an associated, repeatabXe trajectory in the hodograph 
pXane there ought to be some deviation from pureXy 
stochastic behavior. • 

It v 70 uXd be usefuX to formuiate quantities sensitive to 
such subprocesses. There are two cumuXatxve quantities that 
come to mind; the integraXs 

' CO ■ , , , 

/ / £-< dq> ? (u,v) dudv ( Vi a) 

and 

P{u,v) dudv (lib) 

The magnitude of the iatter integraX xn comparison with the 
former provides a quantitative measurs of how much swirX 
is present among the drift vectors. 

In addition to being possibXe pattern recognition 
tooXs, integraXs (11) have an interesting physical 
interpretation in terms of a qeneraXized autocorreXation. 
If compXex velocity W(t) is defxned 

VJ(t) = u(t) - jv(t) (12) 

an the compXex, scaXar autocorreXation < m-T'' > (x ) can be 
defined 


: <WW*>(t) = {< UU>“K VV> } (t) + j {< uv>-^ VU> } (t ) (13) 

I ’■■■ • • 

The one dimensional autocorrelation <uu> (x) can be expressed 
directly in terms of measured forward differences <du>{U/V) : 

i , ■ , ' ■ , , , ■ 

* ■ ■ ■ ■■ CO ■ ' 

•; <UU> (t) = // u {u-Kdu> (t) > P(u,v) dudv (lA) 

-txitis stsp ^ is ths tiiti 0 Qsiay» Fjroin eQJstions (13) 


land (14) 

I ' 00 . 

I < W'J*>(r ) - // g.-t< d 2 > (t) }P (u/v) dudv 

I —00 

i , . ' ■ ■ ■ ■ ■ ■ ■ ■ 

j + j //c[x{ a-K dg^> (t) } p(u,v) dudv 


= ( (u')“ + (v’)^+ //q*< dg^> (x) P(u,v) dudv) 

—CO ■ 

CO 

I -r j //qsK dc[> (x) P (u,v) dudv -(15) 

j : ■ ■ 

■Equation (13) is related to the ordinary two dimensional 
'time correlation tensor 

I . ■ - : • 

1- . _ ^ ^ 

I • • R (t ) = < (t) ( 1 6) 

ISpeGifically, 

I ^ < "WT ' > (x ) = Pkk + jC emnS^ri ^ ^ ^ 

Com^patison of equations ( 1 1 ) , (15) / and ( 1 7 > shov;s that tho 

'amount of " sv?irl *' present in tte drift field is directly 
related to the amount of asymmetry in the matrix. 
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• The complex correlation < (t) is a scalar since ii- 

is invariant under coordinate rotation, and it has a 
magnitude and phase, which will be represented S(t) and 5{x) 
respectively : / 


< (t) = S (T)e 


(t) 


(IS) 


I We make the following hypotheses about S and : 

I , ' ' ' , , 

I ' 

j S(t) is a measure of the correlation between 

i : ^ ^ 

i coherent structures, and 

j ■ • 

I b. ^ (■t) is a measure of the local rotational 

i distortion of the velocity field. 

The physical basis for these statements is illusurated it. 
Ifigure 9. The figure shows the complex correlation 
[(normalized by < u2>-k v^> , the mechanical energy) measured at 
|the location x/D = .3 and r/D = .5. This is the position in 
ithe jet where the initially laminar shear layer has rolled 
!up into regularly spaced, discrete lumps of vorticity v/hich 
(give rise to periodicity in the u and v autocorrelations.^ 
The data points shown in the figure are evaluated at 
multiples of the digitization period. The correlation i^ 
real and unity at zero time delay, and traces out a decaying 
spiral .v;ith increasing t . Pnasa angle * (t ) is a 
monotonically increasing function of time delay, which 
presumably reflects the large local shear, while magnitude 
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S(T) is a generally (although not monotonically) decreasing 

knction of t . It should be emphasised that there are t^« 
independent time scales present: the time scale of the 

'periodicity, which is the time for the spiral to malce one 
Uvolution, and the time scale associated with the decay of 

;s(T). This latter scale depends on how well the discrete 

Slumps of vorticity are correlated with each other. 

,l ■ •' ' ’ ■ ■ ■ 

I The limitations of the concept should be examined. 

'There are two extremes where the complex correlation 

Lupplies little new information over ordinary u and v 

Fir-sfc, if the flow is periodic and 

jcorrelated over many periods, as in figure 9, then the 

Lutocorrelations appear as damped cosines, and their 
'envelope contains essentially the same Information present 
L the decay of S(x) . Second, if the flow is well-developed 
Iturbulent flow, the u autocorrelation is always positive and 
contains most of the information present in the S(t) decay. 
Between these extremes, where the autocorrelations are 
periodic but highly damped, the complex correlation is a 
useful, and consistent method of measuring the correlation 
over the field as a whole, rather than over a single 

structure. 

The measured forward differences are evaluated at one 

aelay time, namely the 
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useful to be able to use that vector field to estimate the 

i 

for\>?ard differences for longer time delays in order to 
estimate the correlation at arbitrary Whatever scheme is 

devised/ there are two constraints that serve as confidence! 

j ■ 

I ■ * ■ ■ ' , » 

(Checks on the estimate. First the magnitude S(x) of the 
correlation must generally decrease (usually monotonically) / 

and second/ the time averaged forward difference must be 

i .■; , ■ ; 

zero / 1 

t ■ '.y 

5 ■ ' ' ' ' ' 

I //<d£>(u/V) P(U/V) dudv = 0 ( 1 9) 

I ' ^ ■ • ' ■ 

■in order that the mean (< u> /< v> ) be time invariant. ; 

!■ : ■ ■ . ■ . '■ ; 

' [ 






i 
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t 

i ^ ^ . 

I ^ 

■ 4. PRELIMINARY JET SURVEY 

i ■ 

t * 

.'4 . 1 Mean velocity profiles 

i ■ 

: In order to establish that the model jet is similar to 

those of other experimenters/ radial profiles of mean 
Velocity were measured with a hot wire probe and traverse at 
Wxial positions of 2/ 4/ 8/ and 1 6 diameters downstream from 

;the nozzle exit plane. The m.ean velocity profiles/ along 

i 

yith profiles of rms turbulence are shown in figure 10. 

I.:' 

I decay with axial distance Oi. centerline mean 

ivelocity is plotted logarithmically in figure 11. Note that 
fhe asymptotic 1/x decay begins about eight diameters 
Idownstream. Data from Crow and Champagne (1 971 ) exhibits 
good agreement. 

1 Linear shear layer spreading / observed along the 

^potential core region of the jet/ is evident in figure 12. 
(The local shear layer thickness h (x) / derived from the_ mean 

‘profiles/ is defined by the boundaries <u>/Do = .10 a 

f u> /Uo = - 9 9 - Note that the initial sheaf layer t^ 
at the nozzle exit is .05 diameters. i 

:t ' ' " - - ' ■ ^ : ' - ' ' - \ '■ 

The self~preservation of mean velocity profiles is 
examined in figure 1 3 . The similarity parameter n 4-S 



F i g u r e 10. P ro f i 1 es of mean ve 1 oc I ty = - t u rbu 1 en t 
velocity. Closed symbols rsfer to mean 
veloci ty. and open symbols refer to 
turbulent velocity. , 
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defined: 


■ ^ = _£Z2. - jlL- (20) 

i x/D 

Reasonably good collapse is observed as far downstream as 
eight diameters. - 

.4 .2 Mean turbulence profiles 

; The profiles of rms fluctuating axial velocity, plotted 
on figure ■ 10, were measured v/ith the hot wire probe at the 
!same time as the mean velocities. They are compared with 

I ... ■ . . 

turbulence profiles measured by other investigators, and 
>'/ill serve as a comparison for the joint velocity 
measurements presented in chapter 6. 

j The turbulence level measured at the n,oz 2 le exit plane 

t - ■ ■ ■ . 

'is shown in figure 14 over a range of exit velocities. For 
most of our measurements, an exit velocity of = 1 00 f os 

v/as seleGted; this corresponds to an exit plane turbulence 
level of .3^. 

Self-preservation of turbulence profiles, figure 15, is 
evident at least as far downstream as eight diameters. 












5. NEAR FIELD MEASUREMENTS 


5.1 Near field static pressure 

There is alv/ays soma question as to V7hat a microphone 
measures in the presence of mean flow. Work by Siddon 
(1 969) ■ and later by Fuchs (1972) indicate that there are 
probe configurations that allov; microphones to respond 
mainly to static pressure fluctuations. For the. 
measurements to follow the ‘ microphones v/ere placed just 
outside the edge of. the jet v;here mean velocities V7ere 

ismall, and v/ere oriented so that the plane of the diaphragm 

* ■ 

was roughly parallel to the mean flow. 

i - 

i ■ 

I Near field static pressure fluctuations can be 

f ■ ■ ■ ‘ ■■ ■ 

•expressed as a voliame integral over momentum fluxes inside 
;the mixing layer. The consequent spatial averaging makes 
■static pressure particularly sensitive to large scale 

f . . V . . . ", ' 

■coherent structures, as indicated by Lau, Fisher, and Fuchs 
; (1972).. t- ^ ^ 

i ■ . ' , - ■ - 

5.2 Near field time scales 

i In order to characterize the mean shear layer 

structures from their near field pressure signatures, 
varidus .statistical methods v/ere used to derive 


characteristic time scales. It was possible to measure mean 
convection times from crosscorrelations between two 
microphones, v/hile information about the time scales of the 
turbulent structures was derived from autocorrelations, and 
from wave counting. The microphones v/ere located on the 
surface of an imaginary cone with a 10° half -angle and 
intersecting the inside diameter of the nozzle exit. Auto- 
and crosscorrelations of microphone signals ware measured 
using the digital correlator, v/hile a counter and 
oscilloscope V7are used for the wave counting. 

The measured crosscorrelations are sho\-?n in figure 16 
for various downstream positions. The two microphones were 
separated 5.0 mm axially; the time delay tc corresponding to 
the crosscorrelation peak is taken, to be the average 
convection tim.e between the microphones. The mean 
convection speeds, defined as separation distance divided by 
convection time, appear in figure 17 as a function of 
downstream position. The convection speeds are normalized 
by the; jet exit speed,: and the ratios are in good agreeraent 
with previous msasurements along the potential core- 
Downstreara of the core the Gonvection spaed decays 
asymptotically as 1/x v/ith a ntagnitude of about half the 
local centerline speed . 

The autocorrelations measured at various 


dovmstream 
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positions are shown i„ fi3^re 1B, nomalizeS by the rn.s 
-gnax. interpretation Of tha associated tin.e scales is 
n,ore difficult than with the crosscorrelations. Along the 
potential core (x/D < 4 ) the autocorrelation appears like a 
damped cosine, a.nd the associated time scale should be th. 
period Of the mean structure. Por a couple of diameters 
downstream of the core there is no characteristic time scale 
at all, while even farther domstream Cx/d > 7, there 

reappears a region ^ of anticorrelation. Since there is no 
periodicity in these latter autocorT-elatiors 

'-he associated 

txme scale is a measure of the coherence of the mixing layer 
turbulence at that point. The axial distribution of the 

time scales as they have been defined appears in figure 19 , 

normalized by exit speed and jet diameter to fora an inverse 
Strouhal nurher. x„ 

.multiplied by the local convection speed to provide as 

estxmate of the coherence length. The _ 

scale appears 

to grow linearly with downstream distance i- }-eeo1nu , • .w 

, ■ — - -f^eepxng v/xth 

observed linear mixing layer spreading. 


^ "cxiti^ soc, Tnoxs 

closeXv 

J between successive 

coherent 

- mixing region of the 

D / 3. 

■ -- ..ficrophone was 

iTs^iOii/. Hixct 

triggers 


used _ to produce, a narrow pulse ea 
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Figure 1 8 . Near field a u toco rre 1 a t Ton s . 
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signal passed through zero volts v/ith a positive slope. The 
pulses V 7 ere counted on a digital counter to derive an 
average passage frequency. In order to reduce the influence 
of the fine scale turbulence the microphone- signal was low 
pass filtered an octave above the measured frequency and 
when the wave counting v/as repeated there v;as negligible 
change in the results. The measured passage frequencies/ 
normalized as Strouhal numbers/ are shovvn in figure 21 . 
Figure 22 shows the local con'/ection speed divided by the 
passage frequency to form the raean spacing X. The spacing 
appears to increase linearly along the potential core but 
farther downstream (x/D > 7) the spacing asymptotically 
approaches a constant of about three diameters . 

Comparison of figures 20 and 22 implies a 
contradiction: the size of the structures eventually becomes 
larger than the spacing betv/een them. The statistical 
implication is that the pressure signal is increasingly 
skewed towards smaller time scales v/ith • inereasing 
downstream distance. That is, the most probable time scale 
(reflected in the autocorrelations} is increasingly larger 
than the average time scale. It is possible that this 
reflects the pairing process in that interactions between 
adjacent concentrations of vorticity become increasingly 
intermittent downstream of the core. On the other hand the 
one-to-one correspondence betvresn lumps of vorticity and 
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near field pressure pulses may simply disappear v/ith 
dov7nstream distance. That interpretation v/ould be 
consistent with the azimuthal correlation measurements shown 
in figure 23. The outputs of tv70 microphones placed on 
either side of the jet at the same dov^nstrean distance were 
crosscorrelated, and the value of the maximum 
crosscorrelation coefficient is plotted against downstream 
distance. The correlation is negligible by the end of the 
core, in qualitative agreement v/ith Mollo-Christensen ' s 
measurements (1957). Again this can be reconciled with the 
pairing model in that the vortex rings are observed to lose 
their axisymmetry through the production of azimuthal modes. 

i 

i From the considerations above it is apparent that near 
field pressure measurements are open to a lot of 
misinterpretation. Consequently most of this investigation 
is directed towards mixing layer measuremen.ts . 
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6. MIXING LAYER STATISTICS 


6.1 Joint (u, v) statistics 

I ' 

; . For the purposes of this section the jet was considered 

to consist ■ of three separate regions: the region close to 
the nozzle exit (x/D < .5) v/here the shear layer undergoes 

transition from laminar to turbulent spreading; the 
potential core region C.5< x/D < 5) v/here the cylindrical, 

turbulent shear layer gro%'7S linearly until the jet is fully 
mixed; and the transitional mixing region (5 < x/D < 20), 

which has been observed to be the major soxrrce of acoustic 

I • ■ 

radiation. Data w^s collected at representative locations 
in each region as summarized in table 1 . 

I Data from each location v;as digitized and analyzed to 

i ■ ■■ ■ , , , ■ ■ , , , ■ • ■ ■ ■ 

produce joint probability densities, drift vectors, and 

diffusion tensors. The results will be presented in the 

following pages. Joint probability d.ensities are displayed 

\ 

in terms of equally soaced contours of constant probability . 
For clarity the presentation of the drift vectors is 
modulated by the corresponding joint probability density. 
The vectors are normalized so that the maximum vector 
magnitude is one cell width. Diffusion tensors are shov;n 
diagonalized in the coordinate system of their principal 

■ . 62 . / 
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axes. They are similarly modulared and normalized. In each 
case- the origin of the coordinate system is the mean 
velocity (<u>/<v>). Where possible zero velocity is 
indicated on the joint probability figures by a tick mark. 
Four straight lines radiating outwards from the mark depict 
the calibration limits. The inner pair of lines define the 
actual domain of the calibration, while the outer pair 
depict the angles of the hot v/ires. 

Note that the time step is fixed at 40 microseconds for 
all the drift and diffusion measurements. This value of dt, 
which is just the digitization period, corresponds to a 
resolution bandwidth' of 12.5 kHz . The power bandwidth of 
the turbulence is less than 5 kHz. 
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Region 

Laminar/ turbulent transition 
Potential core 

Transitional mixing 


x/D r/D n 


.2 .5 0 

.3 .5 0 

2 .5 0 

.6 .05 

4 .4 -.025 

.6 .025 

8 0 -.062 
.5 0 

T.O .062 

12 .5 0 


Joint statistics: 


Table 1 . 

schedule of measurement locations. 


The following parameters were generated to describe the 
joint statistics at each point; 

a) Rms axial velocity fluctuation 

b) Ratio of axial to radial turbulence u'/v' 

c) Mean shear stress uv - u v /u'v' 

d) DRIFT // /< dq> |/ dq> P(u,v) dudv 

e) DIFFUSION / P(u,v) dudv 

The variations of u'/U^, u'/V/ and {< uv>-c ux v> }/u* v* with 
axial position is summarized in figure 2 4 along the ray n = 
0, while figure 25 shows the variations across the snear 
layer. The measurements of other investiga.uors are included 
for comparison. Although there is considerable variation 
between investigators, our measurements seem to fall v;i thin 
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(a) Profile of u'/v' ratio 



(b) Profile of Reynolds stress 
Ra d T a 1 p r of i 1 e s o f mixing lay e r s t a t i s t i c s 
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the scatter. The axial variations of DRIFT an<i 
DIFFUSION/DRIFT are shown in figure 26 , v/hile the 
corresponding variations across the layer appear in figure 
27. 

Recall that the ratio of DIFFUSION to DRIFT is a 
measure of how "noisy" the signal is; the smaller the ratio 
the more deterministic the flov;. The implication is that 
statistically the turbulence is increasingly noise driven at 
increasingly large downstream positions ^ and at increasingly 
large radial positions . 

6.2 Laminar/turbulent transition region 

The flow leaving the nozzle exit is initially laminar. 
The shear layer that develops at the edge of the jet is 
spatially unstable and the resultant v/aves grow with a 
Strouhal nxiraber of .2 (based on shear layer thickness at the 
nozzle exit) v/hich is in agreement with the "eigenfrequency" 
of the most amplified spatial instability Ce.g. Lessen, 
1949). At some downstream position the waves coalesce and 
roll up into cylindrical rings of concentrated vorticity 
which interact \vith each other. From corroborative data the- 
first pairing is knov/n to occur at an x/D bat%-resn .2 and .3. 
The corresponding joint probability densities > P (U/ v) , are 
shown in figure 28. From figures 24 and 23, it is apparent 
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that the aspect ratio u'/v' of the joint probability 
distribution decreases dramatically after transition. This 
can be understood after examining a short sample of the u,v 
velocity record figure 29. At x/D = . 2 , the signal is 

^ i 

dominated by low frequency "whipping" of the whole jet so 
that the statistics reflect the rather large mean shear at 
that point. At x/D = .3, the signal is dominated by the 
strong vortex structures associated vrith roll up and 
pairing. Consequently, the flow process appears more 

deterministic. ■ 

■ ■ ■* 

Measured drift vectors <da> (u,v) are shown in figure 
30, while figure 31 shows the measured diffusion tensors.; 
The lengths of the double vectors are the eigenvalues of the 
tensors while their orientation represents the axes of the 
corresponding eigenvectors. 

The small DIFFUSION/DRIFT ratios, and large anisotropy 
evident among the diffusion tensors reflect the 
determinateness of the flow at these small downstream 
positions. There is considerable "swirl" associated with 
the drift vectors in the sense that they seem to align vrith 
the probability contours to form closed trajectories in the' 
hodograph plane. Tlais behavior seems to be associated with 
the periodic nature of the flov; field at this point, as 
discussed in section 3.4. 
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6.3 Potential core region 

Measureitisnts were male mio-way along the po L-sntial core 
(x/D = 2) , and towards the end of the core (x/D = 4) . Tv;o 

radial positions were selected for each axial station. 
Measured joint probability densities, saiuple velocity 
records, drift vectors, and diffusion tensors are presented 
as figures 32-35. 

The figures above, in conjunction vxxtn rxgures 24 - 2 7, 
indicate the statistical parameters u'/U and u' A’ are 
nearly constant with dov^nstream position %v^hile the shear 
stress slowly increases to a maximum near x/D = 4. DRIFT 
decreases v;ith dov/nstream distance wnxle the ratxo 
DIFFUSION/DRIFT increases roughly linearly v/xrh dxstance . 
All parameters except DRIFT increase across tne shear layer 
away from the jet axis. However DRIFT is a measure of flow 
acceleration and will necessarily decrease as length scales 
increase. 


For positive values of the similarity parameter, g, the 
joint probability densities (e.g. figure z>2o) exhibit an 
"edge" in the third guadrant characterizec. by parallel, 
closely spaced contours. This effect becomes more 
pronounced v?ith increasingly positive n . The drifts are 
relatively small near these "edges" while the diffusion 
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Figure 32. Concluded . 
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tensors become small and anisotropic, v?ith major axes 
aligned parallel to the edge (e.g. figures 34b and 35b) . 
The significance of the edges is unclear, and may simply be 
an effect of the calibration. The edges lie along angles in 
the hodograph plane that are generally outside of the 
calibration domain. 

6.4 Transitional mixing region 

Since the most intense noise producing region of a 
subsonic jet has been measured to be aroond eight diameters 
downstream from the nozzle, this region v?as examined in 
detail. Data v;as taken on the jet axis (r/D = 0) , at the 
radial position of maximum turbulence (r/D = . 5), and at the 
outer edge • of the jet (r/D = 1.0), where the turbulence is 
intermittent. For completeness additionaL data was recorded 
at (x/D = 12, r/D - .5). The joint statis tics as well as 
velocity records are shown in figures 36—39. 

Statisticallj? there is nothing strikingly unique about 
this region except that the mixing occupies a greater 
volume . For the most part the statistical parameters are 
similar in magnitude to those in the potential core region, 
and decay slov/ly with dpv/nstr earn distance (xigures 24 and 
26) . The exception is the ; DIFFUSIOJi/DillFT ratio Miich^ 
continues to increase linearly with x/D . Consequently at 
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some point In the transitional nixing region the velocity 

signal can be modelled statistically as a noise driven 
process. . . 

6.5 Generation of complex correlation 

The complex scalar correlation < ww*> (x) was defined and 
discussed in section 3.4. The correlation was then measured 
at position (x/D=.3, r/D=.5) according to the scheme 

described. The periodicity of the flow at that location is 

reflected in the measured correlation, figure 9. it. is of 

xnterest to measure the correlation for locations farther 
downstream where the periodiGity has disappeared. 
Accordingly, figure 40 shows the complex correlations for 
positions (X/D = 2,8; r/D = . 5) . These correlations behave 
quite differently from the earlier one. The imaginary part 
of the correlation is small, and the change in phase angle 
is no longer monotonic. This behavior in phase angle is 
difficult to interpret. 

The drift vectors ar^ ‘^ifferenGes 

orientations are maasures of vmiGh direction the fluid was 
traveling in space during the time increment 
counterclockwise solid body rotation in a spatial coordinate 
system being swept past the probe would appear as clockwise 
trajectories in the hodograph plane. 




direction of change of phase angle observed in figure 9 is 
consistent with the local shear. In figure 40, however, the 
phase angle changes in the opposite direction. Possibly the 
local fluid distortion is dominated by strong intermittent 
influences such as .entrainment rather than by mean shear. 

The decay of the correlation magnitude S(t) at two and 
eight diameters is now clearly monotonic in bime delay. The 
decay of S(t) is presented in figure 41 as a function of t 
normalized like an inverse Strouhal number. All three 
locations (x/D— .3,2,8; r/D=.5) are included and the display 
is essentially log log S(t) vs log ( t) . Any exponential 
behavior of the decay will appear linear under that scaling. 
At the location (x/D=.3, r/D=.5) there are oscillations 
superimposed on -the decay. However the dafa seems generally 
to fall off like exp (-/t) . At the other positions, however, 
the decay was very closely Gaussian, for small , and 
generally approached something like a simple exponential at 
larger delays . 

From the observed decay of S (x ) it is possible to 
define a correlation length scale . If t ^ is chosen as a 
characteristic time scale such that S (x c ) /S (0 ) = 0 . 1 , then. 
X cUc/h is the convected correlation length scale normalized 
by the local shear layer thickness. This length , scale 
appears in figure 42 for various downstream locat.lons . For 
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the position (x/D=8, r/D=.5) the local jet radius v;as chosen 
as the shear layer thickness. The correlation length scale 
decreases from about T1 shear layer thicknesses at (x/D=.3, 
r/D=.5) to about 0.7 shear layer thicknesses at subsequent 
axial positions. 

It is an interesting question whether there is enough 
information in the drift field measured over one time step 
to estimate the drift field and hence the complex 
correlation < W*> (t) for arbitrary time delay x=nTQ . 

TThat ever scheme is selected, there are constraints on the 
drift field as well as on the correlation itself that serve 
as checks on the reliability of the estimate: 

1 ) as the time step t is increased the drift vectors 
must increasingly point tov/ards the mean {< u> ,< v> } , 

2) since the statistics are stationary, the integral 
//< dq> (u,v) P (u,v) dudv must add up to zero, 

■ 3) from physical arguments the magnitude S (t) must 
generally decrease. 

The simplest scheme that comes to mind is to associate 
a direction and magnitude with each cell (u,v) in the 
hodograph plane based on the measured drifts, to note in 
v/hich cell the tip of each < do> vector falls, and to then 
add the corresponding magnitude and direction. According to 
this scheme of course the time delays are restricted to 


The vector 


integer multiples of the original sampling rate, 
arithmetic is shoim schematically in figure 43. The 
resulting estimate of the complex correlation is shown in 
figure 44 for the position. (x/D= 8, r/D=. 5) . The correct, 
measured correlation appears in figure 40b. Comparison 
between the tv 70 indicates that the scheme doesn't provide a 
good estimate. Unless- a more elaborate estimation technique 
would provide more accurate results, the most .efficient 
method of measuring complex correlations is to use a Fast 
Fourier Transform algorithm and measure the individual 
elements of the tvvo—dimensional correlation m.a-crix . 








7 . EVIDENCE FOR PAIRIJ^G HYPOTHESIS 

7.1 Conceptual basis • 

The role that vortex pairing plays in the growth of a 
turbulent shear layer has been demonstrated by Winant and 
Browand (1974) . Vortex pairing is considered to play a 
similar function in a turbulent circular jet. However the 
mechanism is complicated by the fact that unlike the free 
shear layer the jet has two length scales: the exit diameter 
as well as the local shear layer thickness. Certainly in 
the early stages of growth the jet shear layer is similar to 
a free shear layer. In fact since a cylindrical vortex 
sheet lacks a most amplified mode of instability, the 
laminar instability is probably tv/o-dimensional . Analogous 
to a two dimensional shear layer, the jet instability rolls 
up into Ixrnps of concentrated vorticity which coalesce into 
rings of vorticity which interact to cause the initial 
mixing layer spreading. This much is evident from the flow 
visualization studies {e.g. Browand, Ghu, and Laufer, 1 975 ) 
at least for low and m.oderate Reynolds numbers. The nature 
of tne interaction between turbulent vortex rings is 
suggested by the kinematiGs of vortex rings idealized from 
potential flow theory and is indicated diagrammatically in 
figure 45. Consider the case of a ring where all of the 
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vorticity is confined to an annular core of radius a which 
is small in comparison v;ith the ring radius R. Then the 
speed of the ring depends only on R and the circulation r 
(Batchelor, 1967): 

U_ = r ln(R/a) (21) 

Wr • . 

Note that for constant circulation the ring moves faster as 
the radius becomes smaller. . Now consider tv;o rings one 
behind the other: the induced flow from the forward ring is 
radially inwards tending to compress the rear ring, causing' 
it to accelerate, while the flow induced by the - rear ring 

causes the forv/ard ring to spread out and decelerate. The 

» 

rear ring will overtake and pass through the forv/ard ring. 

According to potential flow theory this ’'threading” process 

would repeat itself indefinitely. In practice the cores of 
the rings are observed to coalesce producing a single ring 
with twice the circulation. Unlike the free shear layer 
there are two scales involved in ring pairing: the ring 
diameter Ci.e. jet diameter) and the core diameter (i-e. 
shear layer thickness) . It is an open question whether 
pairing continues v;hen the core has grown to a large 
fraction of the jet radius . The flow visualization studies 
do not help here because there is so much mixing that the 
dye or hydrogen bubbles are ooscured. 


In the sections to follow, eduction experiments are 
described that are designed to detect pairing, and to detect 
its signature in the near field. It is an easy matter to 
detect the roll up of the laminar instability and the 
occurrence of the first pairing. It is more difj.icult ro 
detect subsequent pairings, although our attempts are 
described. Finally time averaged near field measurem.ents 
are interpreted according to the pairing hypothesis to model 
the mean location and variance of the pairings. 

7.2 Correlation maasurements inside shear layer 

The laminar shear layer instability is suj-j-icientiy 
airiEilified at about x/D = 0.2 that it can be easily detected 
V 7 ith a hot wire probe. The probe was placed into the shear 
layer at the radial position corresponding to maximum 
turbulence level and autocorrelations of the sa.gnal over a 
range of downstream positions are shown in figure 46. For 
axial stations less that 0.4 diameters the autocorrelations 
appear like damped oscillations with constant frequency . 
Normalized by the initial shear layer thickness of .05 
diameters (from figure 12) , and velocrry u = Ug/2, tnis 
frequency corresponds to a Strouhal nimber or 0.2. 
Two-dimensional linear stability theory (e.g. Lessen, 19'-i9) 
predicts the Strouhal number of the pre ^ rentxally amplifxed 

- • ' . ■ ■ ■ ■ ' \ . /r . .. . ■ ' . ’ ■ ' ■ ■;/' - 

mode to be between .2 and .3. So the data is so far 
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Figure kS. Tu rbu 1 en ce a u tocorrel a t ions inside 
mixing layer. fr/D = .5) 
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consistent with the model. 

There is sufficient temporal coherence at these early 
upstream positions that there is evidence of the first 
pairing in the autocorrelations. Betv/een stations x/D=.4 
and x/D=.7 every other v?ave in the autocorrelations 
disappears. This amounts to a discrete shift in frequency 
to the exact subharmonic of the eigenmode. The implication 
is that there is a discontinuous doubling in length scale 
somewhere between x/D=.4 and x/D=.7. This is certainly 
consistent v/ith vortex pairing. Beyond x/D = .7 the 
autocorrelations exhibit a continuous increase in the time 
scale associated with the periodicity. it turns out 
there is at least one more pairing that can be distinguished 
using conditional sampling. The continuous shift in 
frequency observed in the autocorrelations is the result of 
the statistical variation of the location of the subsequent 
pairing (or pairings) . 

7.3 Two point eduction techniques 

In order to derive some spatial information about the 
first and subsequent pairings, simple two point eductions 
were performed across the shear layer at various axial 
locations as far downstream as x/D=9 . Two constant 
temperature single wire probss were usedr the triggering 
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probe was placed in the core next to the shear layer and the 
probe whose signal was being educted vras placed in the shear 
layer at the same axial station as the triggering probe but 
at varying radial positions. A pulse was generated vmenever 

the trigger probe signal exhibited a positive slope and a 

level above the mean greater than one r.m.s. of the 
fluctuation. This condition is a standard first attempt in. 
conditional sampling of turbulent shear flows (e.g. 
WillTnairth and LU/ 1972) . 

The educted waveforms appear in figure 47. Each 
subfigure corresponds to a particular downstream station , ^ 
and in each case the indicated radial positions are 
relative, with positive r pointing in the direction of the 
entrainment region. Closest to the no.zle (x/D=.20, figure; 
47a) there is an abrupt 180° phase shift in the center of 
the shear, layer accompanied by small fluctuation amplitudes . 
This behavior is suggastive of. a simple vortex sheet. 

Another 0.1 diameters downstream; however, the behavior has 
changed. The phase shift is spread out over 15 of the 

shear layer thickness, and in the center there are twice as 
many crests (figure 47b) . The corresponding harmonic 
frequency is absent in the power spectra (e.g. figure 46) , 
which indicates that instead of a frequency doubling there 
is a double valued phase relationship between the trigger 
condition and the educted waveform at this radial location . 
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Figure 47. Educted veloci ty signatures across 
mixing layer. 
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(b) x/D = .28, Ar/h = -.13 to .07 


Fi gure A/ . Conti nued . 
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(d) x/D ^ 4.0, Ar/h = -.23 to .46 
Figure 47 • Continued. 
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This double valued phase relationship is observed in, the 
eduction measurements of Winant and Browand. (1974) in the 
free shear layer and is caused by the pairings. The double 
structure is barely discernable at x/D=.9 (figure 47c) v;here 
the associated time scale is the second subharmonic of the 
eigenmode. The double structure here probably corresponds 
to a third pairing. At farther dovnstreaia positions, x/D = 
4 and 9 v;hich are positions of interest acoustically, it is 
unclear whether the waveforms represent a phase shift across 
a coherent structure or an increasing phase delay due to the 
variation of mean velocity across the layer. 

A more sophisticated trigger condition must be 
developed before pairings can be detected farther dovmstr earn 
than about a diamater. Since the ultimate intention is to 
be able to detect the pairings from the near field, a 
reasonable first step is to generate -various kinds of 
conditionally averaged hear field pressure waveforms based 
on different sampling conditions with the goal of selecting 
a subprocess suggestive of a vortex pairing . Accordinglv , 
an x-wire probe was placed in the shear layer at the 
position (x/D = 2 . 0 , r/D=. 60) , which is slightly tov/ards the 
entrainment side of the turbulence peak and at an axial 
station where -the Goherence across the get is still hieh . 
Two microphones v/ere placed in the near field of the jet at 
the same downstream location, x/D=2.0. One was positioned 
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immediately next to the x-wire probe and the other 180° 
opposite. The turbulence and near' fields v?ere recorded 
simultaneously, digitized, and various sampling conditions 
derived from the u,v record were investigated. On the 
whole, the attempt was not very successful. The main 
difficulty is that the jet, unlike a boundary layer, is not 
constrained by fixed surfaces, and the spatial wandering of 
the jet probably obscures the subprocesses of interest. The 
types of trigger conditions investigated and the results are 
siimmarized: 

1) Various kinds of VITh (Variable Interval Time 
Averages) averages of the sort pioneered by Kaplan and 
Laufer (1969) in connection vrith the turbulent sublayer 
were applied to both the u signal and the uv product 
signal, but produced no interesting educted signatures. 

2) Positions (u,v) v;sre selected in the hodograph 
. plane that were associated with trajectories of high 

turbulent kinetic energy, and educed signatures were 
computed for the uv record as v/ell as for the pressure 
records . The trajectories v;ere computed by looking 
forward 1 0 time steps and backwards 10 time steps 
whenever the u,v record occupied the cell (u,v) being 
investigated. The kinetic energy was averaged 
^ v/ith the trajectories . IvO consistent 

trend was observed in either the educted pressure 
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signatures or in the educted uv signature . 

31 Educted u and v signatures along V7ith educted 
pressure signatures are shown in figure 48, and were 
educted firon the triggering condition 
(uv - <u><v>)/u’v' greater than 1.5. The upper two 
curves are the educted u and v signals, the center 
curve is the. ensemble averaged uv product, and the 
lower tv70 curves are the near and far microphone 
signals. Note that there is a clear phase relationship 
between the conditionally averaged pressure signals and 
the condition of large uv product. However the time 
scale of the educted pressure waveform is clearly that 
of the mean structure, which implies that the trigger 
condition is sensitive to the mean structure rather 
than interactions between structures. This outcome is 
typical for the types of triggering conditions 
investigated . 


7.4 Estimation of mean pairing locations 

As things now stand the pairing process, clearly 
observable by visual technigues, can be detected only in the 
earliest stages of shear layer spreading. More 
sophisticated sampling techniques such as u, v space 
statistics have been ineffective in identifying the 
occurrence of pairings farther downstream. 
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Figure 48 
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It is instructive nevertheless to model 
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the near field 
measurements according to the pairing hypothesis to deduce 
V7hat the mean locations and variances of the pairings \-7ould 
have to be. If the near field time scales represent the 
periodicity of vortex ring passages, then the mean locations 
of the pairings can be deduced deduced directly from the 
time scales using conservation of waves. Within a given 
interval the rate, h, at V7hich vortex rings 
disappear is 

h = f Cx^) - f(x2) (22) 

V7here f (x) is the frequency of vortex ring passage past 
position X. If X. is the mean location of the '’i"th 
pairing, then half of the rings passing x=X . will have 
disappeared due to pairincs before reaching x=X.^^ , the mean 
location of the "i+1’’st pairing: 







( 23 ) 

From the 

two-point 

eductions 

across the 

shear 

layer it 

appears 

that the 

first pairing takes 

place 

at about 

x/D = .2. 

Then based on the 

V7ave counting 

measurements, the 


mean vortex passage frequency as a function of x is: 

f(x) = f^/x : i2h} 

It follows from equation (23) that 



ORIGINAIl PAGB JB 

OP POOR iftUALira ■ Jhb 

X;^, - 2X. (25) 

So based on this model the mean positions of the first seven 
pairings are suinrriarized in table 2. 

Pairing # Mean location 
n Xn/D 


1 .2 

2 .4 

3 .8 . 

4 1.6 

5 3.2 

6 6.4 

7 12.8 

Table 2 . Distribution of pairing locations 

Since directional microphone measurements by Chu, Laufer, 
and Kaplan (1975) have localized the noise produeing region 
of subsonic jets v;ithin the range four to ten diameters 
downstream, the point of this exercise is that two or at 
most three pairings would have to be responsible for most of 
the radiated noise if pairings are in fact the important 
acoustic sources. 

According to the model, vortex pairing is the mechanism 
for the sprsading of turbulent sheaf layer s . The fact that 
the spreading is observed to be linear permits an estimate, 
of the the variance associated with pairing locations. For 
the purposes of this crude estimate, the nth pairing is 
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The resulting shear layer thickness, h (x) , is; 

CO 

h(x) = h^[l+ l 2 ''~'^Aj^exp{-a-X)^/ 2 {aZ )^} d?) (26a) 

w - Itw II M 

n= 1 

An = (/oexp{-{?-Xn)^/2 V (26b) 

2 

5-n “ vortex spacing ^n^j “ ^ ~ \*araance in the 

pairing location as a fraction of the local spacing i?' . 
Figure 49 shows the shear layer growth coitiputed from 
equation (26) for three different values of the parameter o. 
For values of a less than .5 the spreading is not smooth, 
and for values greater than 1.5 the spreading is no longer 
linear. .On the basis of these sinple arguments it would 
appear that the dispersion in pairing location v?ould have to 
be of the order of the mean distance betv/een pairings, and 
that this distance V70uld grow exponentially v/ith downstream 
distance. This large "phase jitter" in pairing location 
relative to the mean structure would certainly account for 
the difficulty in detection.' 
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• 8 . CONCLUSIONS 


1) Correlation methods are adequate for detecting 
the regularly spaced vortex rings along the potentxal core 
of the jet. Simple two point eduction techniques are 
adequate for detecting the initial one or two pairxngs near 
the nozzle exit. Even more sophisticated eductxon 
techniques were unable to detect subsequent pairings from 
their near field signatures. Consequently the role played 
by vortex pairing in turbulent mixing layer spreading in 
moderate and high Reynolds number jets is still in question, 
as is the contribution of pairing to far field noise. 

2) The near field autocorrelations indicate a 
fundamental change in mixing layer structure at the end of 
the potential core . In particular the decrease xn length 
scale and total disappearanGe of near field txme scaa.e imply 

a breakdovm of the regular structures observed upstreaia. 

3) The transitional mixing region of -the jet, v;hich 

is observed to be the major noise source, is not remarkable 
in terms of overall statistical parameters such as 
turbulence levels, shear stress, or statistical ' noisiness 

4) - In terms of the distributive measures of DRIFT 
and DIFFUSION aS defined in this paper, the jet turbulence 
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appears fairly deterministic. Stochastic models which 
consider the jet to be noise driven v/ould be most valid for 
large axial and radial positions . 

5) The hodograph plane description of turbulence is 
rich in statistical information both direct, and. implied. It 
seems liJce a promising pattern recognition tool, and is 
worthy of further development. 
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APPENDIX A 


X-WIRE PROBE CALIBRATION 

A.1 Kings law aporoxima-tion 

One point becomes obvious from the joint probability 
densities measured in chapter 6: the flow spends a 

considerable fraction of time at rather large flow angles. 
This fact makes the calibration procedure a very important 
part of x-wire measurements. 

King's law of cooling relates the voltage signal, e(t) , 
of a constant temperature hot v;ire probe to the velocity 
component normal to the V7ire : 

. _ — — . 

e - A + b/u sxna + v cos-^. (Al) 

where u and v are the axial and radial velocities, a is the 
angle of the 'wire relative to the x-axis , and A and B are 
constants of calibration. Equation (Al) can be inverted, 
yielding: 

2 ' ■" h 

u = a, + a_e + a_e + a,.v ( A2 ) 

,, , , ' ' 1 ■ ' , 4 ■ ' S ■ ' ' 

where the coefficients a. are functions of A, B, and «. In 

I ; 

the case of an x-wire probe where there are two wires 
equation (A2) assumes the general form: 
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+ a . , V 
1 h 


(A3) 


where a., is the ’j'th coefficient of the 'i'th v;ire. 
• J 


In practice, hot v;ire voltages, e., are recorded over a 
range of known u and v, produced by placing the probe in the 
parallel flow region of the jet near the exit plane and then 
yawing the probe through knov;n angles. Enough calibration 
points are recorded so that the set of equations (A3) are 
overdetermined. The equations are least squares fitted to 
the data to determine the ''best fit" coefficient matrix 

(a ). Once the coefficients are determined, turbulent hot 

I J 

wire signals e^ (t) and e^ (t) can be linearized to produce 
velocity records u(t) and v(t) : 




^2!i ■ "ll 

\r rr Y ^ F 

2 1 

^14 " ^24 


where the polynomial 


j = 1 


= u - a.^v 


(A4) 


(AS) 


A Set of measured calibration polynomials Ft are 
•plotted against e . in figure Al , and compared v;ith the 
calibration data (u-a . , v) . “he trouble here is that the 

I 4 

function F(e) always has a single minimum located at the 
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- T and unless calibrati-Oii 

smallest measurea value of (u a.^v), ana un- 
data points can be supplied where Cu-e . ^v) = 0, the 
magnitude of F (e) at its minimum will be greater than zero. 
Unfortunately (u-a^v)=0 amounts either to zero flow velocrty 
or to flow parallel to one of the wires. Either case is 
outside the range of validity of Kings law. This deficiency 
in the calibration is a serious one when its effect is 
assesed in the hodograph plane. The set of inequalities 

(F )min > 0 effectively displaces the origin from (u,v) 
i 

(0,0) to 


= a,JF,)min-a'^(F^)roln 

~ ®2l. ■ “U 

V = (F^)min-(F )min 
• o 2 J 


{ A6) 


'U 




For the calibration polynomials in figure A1 the origin 
would be displaced to (u^,v^) = dO fps, 0 fps) , which can 
be enough to significantly distort the isoprobability 
contours of the joint probability P(u,v). 


A 


.2 Linearization subject to constraxnt 


The problem of displaced origin can be alleviated by 
constraining F. (e. ) to pass through zero at its minimum, 
instead of minimizing th^ the yariances/^^^^^^ t^ 

quantities are defined 


15 ^ 



N 

= i { (u-a . ijV) 

n = l 


n 


F . (e .) } 
I I n 


X .F . (e 

I I 



(A7) 


where the 'n' &'ummation is over the data, X. are 
.undetermined multipliers, and e.Q is the locatxon of 
(F. )min. Each is minimized with respect to the 

coefficients a.j, resulting in a set of Ixnear equatxons 
which vrhen inverted yield solutions of the general form 
(summation over repeated indices not implied) . 




j 


(A8) 


The multipliers X. are determined by substitution into the 
constraint F. (e. q ) min = 0^ In ' practxce xteratxon xs 
necessary to finally determine the set (a.^. ). 

The calibration polynomial obtained in this way is 
plotted in figure A2 along with the calibration data. The 
polynomial is now vzell behaved in that it does pass through 
zero at its minimum. However the fit itself has become 
unacceptably poor, exhibiting a mean deviation of about 4 
fps.. 


A. 3 Linearization subject to correction terms 

The fit can be improved while avoiding the problem 0 . 1 . 

origin shift by adding a correction terra to the Fj (e; ) 

polynomials / equation (A5) . The Gorrection terra adds an 





( a ) C h ;4 n n c 1 1 

Figure A2, Calibration polynomials subject to constraint: 
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extra degree of freedom, allowing the polynomial to pass 
through F=0 without constraint. A clue to tne j-orra of the 
correction term is provided by King's lavr. King's law, 
equation (Al), can be generalized to: 

• U f(6) = e^-Aj (A9) 

where U and 0 are the magnitude and direction of the local 
velocity , and the exponent r is approximately 2 . 
Empirically e^ is always larger than A, so than the serxes 
expansion for equation {A9) is: 

U f (3) = fe"] ' ! (Mgll!: 

B ^ n=0 n ! ( r~n)/l 

consequently , the correction term should be of the order 
e"^. The new calibration polynomial Hj(e.) is defined: 




= u - a . ^V 

I 5 


(Al 1) 


The measured calibration polynomials H. (e .) are shown xn. 

figure A3 along with the data used in the curve fittxng. 

This method produced an acceptable xit (mean c~- .>_<=. u_on 

fps) , and eventually passes through zero. -he slOf- 

becomes quite shallow near H. (e.) = 0. The slope xs steeper 

if a correction term of the order e^ is selected instead of 
e" 2 . Although such a correction term is non motivated by 
King's law it has the same effect as the a “ correctron. 





namely it adds an extra degree of freedom to the curve fit. 
The calibration polynomials are redefined: 


H. (e. ) 

I I 




a e 
13 


a , e 

1 k 


(A1 2) 


= u - a . 

1 p 

The resulting fit shown in figure A4 is the polynomial that 
was finally used in linearizing the x-wire signals to 
produce, the joint statistics discussed in the body of this 


paper . 








